Crossover from adiabatic to antiadiabatic quantum pumping with dissipation 
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Quantum pumping, in its different forms, is attracting attention from different fields, from funda- 
mental quantum mechanics, to nanotechnology, to superconductivity. We investigate the crossover 
of quantum pumping from the adiabatic to the anti-adiabatic regime in the presence of dissipation, 
and find general and explicit analytical expressions for the pumped current in a minimal model 
describing a system with the topology of a ring forced by a periodic modulation of frequency w. 
The solution allows following in a transparent way the evolution of pumped DC current from much 
smaller to much larger uj values than the other relevant energy scale, the energy splitting intro- 
duced by the modulation. We find and characterize a temperature-dependent optimal value of the 
frequency for which the pumped current is maximal. 

PACS numbers: 03.65.Yz, 85.35.Be, 03.65. Vf 



A current with a net DC component can be pumped 
in an electronic system without leads and bias voltages, 
through a "peristaltic" modulation of the transmission 
amplitudes and gate voltages jlj. This effect has both 
classical and quantum components @, Q , and occurs for 
unpaired electrons as well as for Cooper pairs 0, 
When the modulation is adiabatic i.e., when the pumping 
period is much longer than the intrinsic time-scale of the 
system, so that transitions between states do not occur, 
it has long been recognized that the charge pumped over 
a period has a geometric nature 0, [3], Q and is in many 
cases quantized. These geometrical aspects survive even 
in the presence of a coupling between the electrons and 
an external phonon bath, despite the obvious source of 
inelastic effects represented by the bath [J, la, [7| . Unsur- 
prisingly, like in classical pumps, the current in this slow, 
adiabatic regime increases proportionally to the driving 
frequency w/27r, as long as Huj is much smaller than all 
intrinsic energy scales of the system. 

The question we address in this Letter is: of what kind, 
and of what magnitude are the deviations from adiabatic 
pumping that will show up in the DC current when the 
pumping frequency grows higher and higher? What is the 
behavior of the DC current as frequency crosses over be- 
yond the adiabatic and into the antiadiabatic (oj — > 00) 
regime? To obtain specific answers, we shall focus on 
the crossover from adiabatic to antiadiabatic quantum 
pumping (or "stirring" @, H, in a system with the 
topology of a ring, in the presence of dissipation. For a 
particular but reasonable choice of coupling to the bath, 
we find that the dissipative model admits a full analyt- 
ical solution for the steady state current valid at arbi- 
trary frequency. Through that solution we can analyze 



and understand the main predicted features of pumping- 
frequency dependence of the DC current. At low fre- 
quencies the pumped current tracks the known adiabatic 
result, namely DC current increases linearly with fre- 
quency, and the pumped charge is as expected geometric 
in nature (albeit not quantized). However, and this is a 
surprising outcome, the pumped DC current turns non- 
monotonic for increasing u, going through a temperature- 
dependent optimal value and then dropping eventually as 
lo^ 1 for u — > 00. Beyond the strict limits of the present 
model, we also surmise that these results are represen- 
tative of a larger class of orbital doublet systems weakly 
coupled to a generic environment. 



Consider the minimal model constituted by a three- 
site ring as in Fig. [TJ each of the identical sites a, 6, c en- 
dowed with a single nondegenerate electronic level of en- 
ergy €i(t), and different sites coupled by hoppings 7ij(t), 
where i,j — a,b,c. Current pumping can be obtained, 
for instance, by letting 7^ = 70 and externally ac- 
tuating a cyclic variation of the three on-site energies 
e. L (t) = — HA cos (uit + 4>i), with 4> a = 0, <f>i, = —2ir/3, 
4> c = +2n/3. In the perturbative limit (HA <C 70) 
the three-site ring can be replaced by a simpler effec- 
tive orbital pseudospin model obtained by removing the 
totally symmetric state |0) = (|a) + \b) + |c)) /y/3 (of en- 
ergy —270 for A = 0, doubly occupied and irrelevant) 
to retain only the two states \x) = — \c))/y/2, and 
\y) = (2 \a) — \b) — |c))/V6, orbitally degenerate in the un- 
perturbed ring A = 0, with energy 70. The single mobile 
electron now occupies the orbital doublet |x), \y), lead- 
ing to a pseudospin- 1/2 problem with the time-dependent 
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FIG. 1: Schematic of the minimal ring model, realized by 
three quantum-dots a, b, c with externally controlled poten- 
tials e a , €b and e c used for pumping. The 7^ are here inter-dot 
hoppings, generally fixed. An alternative realization could be 
a molecular trimer [Tol . ITl| ] where three atomic orbitals |o), 
I b) and |c) have identical energies, while the hopping integrals 
7ij can actuate the pumping by modulating atomic displace- 
ments. In either case the relevant states of the system can be 
mapped to the pseudospin-1/2 model of Eq. JT|). 



Hamiltonian II 



n s (t) 



HA 



cos(u>t)a z + s'm(ujt)<T x } , 



(1) 



where ct^ are Pauli matrices. The current I(t) = (I a b) = 
— «(?e7o(c^c a — c^Cb) (where q c is the elementary charge) in 
the pseudospin representation is given by I(t) — I (<j v ), 
where I = q e Jo/V^- 

While this is the same type of Hamiltonian previ- 
ously used to study adiabatic pumping we now 
find that this problem allows a more general analyt- 
ical solution for arbitrary pumping frequency beyond 
the adiabatic regime, and in the presence of a cou- 
pled bath, so long as the coupling is weak. The exact 
time evolution induced by (fjQ) can be obtained by not- 
ing that H = (hA/2)R y (ujt)a z R- 1 (ujt), where R y (ujt) = 
e ~iujta v /2^ re p resen ts a uniform rotation by an angle wt 
around the y-axis. Performing this time-dependent uni- 
tary transformation and defining \ip(t)) = R y (u>t)\ip(t)) , 
the Schrodinger equation for \ip) is governed by the effec- 
tive Hamiltonian 



n eS = R-'HsRy-iRyRy = - (Act 2 - Loa y 



This now represents a time-independent field pointing 
in the direction n = (0, — v/v', A/a/), where v' — 
Vu! 2 + A 2 is the associated Larmor frequency. The prob- 
lem thus has a simple solution in this reference frame: the 
spin state \ip) precesses around n, while the current re- 
tains the form I(t) = Io(tfj(t)\a v \^(t)} . The current car- 
ried by the eigenstates |n; ±) ofH e s is Io(n;±\a y \h; ±) = 
^FIquj/lu', respectively. In the absence of coupling to the 
bath, all time dependence of the current is determined 
just by the initial conditions. In particular, the two eigen- 
vectors of 'Hcfr carry a pure (and opposite) DC current, 



while any other initial condition yields a DC plus an AC 
current. The DC component is determined by the pro- 
jection of the pseudospin onto the eigenstates of H c s : 



I = IaP- 



(2) 



with the pseudospin polarization P = — Tr(n • <r ps) ex- 
pressed in terms of the density matrix operator ps in the 
rotating frame. 

Even if one prepares the initial density matrix in a 
pure state, the slightest dissipation will eventually drive 
the system to a different (generally periodically time- 
dependent) steady state. To describe the effect of dis- 
sipation, we introduce the environment in the standard 
form [13| of a heat bath of harmonic oscillators at tem- 
perature T linearly coupled to the charge fluctuations, 
embodied in this system by the two operators ct z and <j x . 
The dissipative part of the Hamiltonian is thus Hb+Hsb 
where: 



= EE 

S,=z,x v 



2m 2 



(3) 



ksb = EEv^F A ^^- (4) 

Here uj v are the oscillator frequencies and A^,, are cou- 
pling constants, for which we assume ohmic spectral den- 
sities [3 Hi J 6 (u) = J2u XlA<*>-<»») = h 2 a^e- u / u -. 

Assuming the coupling to be weak, « 1, and retain- 
ing the lowest-order in ot£, the evolution of the system's 
reduced density matrix ps in the rotating frame is given 
by the master equation [15] 



dps(t) 



-i[H cS , p s (t)} 



-y 



dr 



(5) 



{g c (t) [a*®, Ui{-T)5*{t - T)U (-r)ps(t) 
+ G|(r) [psWk-T)3*(t ~ r)Uo(-r), 5*{tj\ } , 

where a^(t) = Ry 1 ^)^ R y (ujt), and U (t) = 
exp(— iH e ffr/K). The function G$(t) is expressed in 
terms of the spectral density as 



Q(t) 



dtoJ^(ui) 



. . huj/3 . 
cos(cjt) coth — 1 sm(cjT) 



where f3 = (fceT) -1 , and T is the bath temperature. 

When the bath coupling to a x and a z have the same 
spectral density (J x (uj) — J z (lu) = J{oS), even if non- 
ohmic) the form of Eq. ([5J becomes particularly simple, 
since all explicit time dependence disappears and we are 
left with a constant-coefficients inhomogeneous linear dif- 
ferential equation 12j. For — > we find that the sta- 
tionary density matrix ps is diagonal in the basis |n±), 
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FIG. 2: Steady-state DC circulating current / as a function of 
the pumping frequency ui. Solid line: T = current, obtained 
by Eq. ([2} with P = 1. Note that this coincides with the 
purely quantum result (no dissipation) when \ip) — n; — ) is 
chosen for the initial state. Dashed and dot-dashed curves: 
pumped DC current for several temperatures, as obtained by 
the exact expression ((6| for P. 

with a polarization given by 

= (uj'-uj) 2 J + + (uj> + uj) 2 J- 

(uj' - w) 2 c+ J+ + (a/ + w) 2 c_ J_ ' 1 ' 

Here J± = J(ui' ± w) and c± = coth[fi(w' ± cj)/(2fc B T)]. 
The resulting DC circulating current, Eq. ([2]), is shown 
in Fig.[5]for a broad range of frequency and temperature. 

Several comments are in order, (i) In the T — >• 
limit, irrespective of w/A and of the form of the spec- 
tral density, the stationary master equation operator is 
a projector onto the ground state |n— ) of H c s- (ii) For 
uj <C A, we have P — tanh[fi,A/(2fcB7 1 )] as is appro- 
priate for a static Hamiltonian in thermal equilibrium. 
At T = 0, P = 1 and the charge pumped in a period, 
Q p = 2irl/cu = 2irIo/(hA) coincides exactly with the 
Berry-phase result [l6[ of Eq. (21) in Ref. [6J. Never- 
theless, the pumped charge, although geometric, is not 
quantized, and depends on the ratio 70/ A. It can be 
made quite large (though not arbitrarily large) by re- 
ducing the amplitude of the perturbing field A in the 
range frus <C HA <C 70 ■ (iii) Finally, in the antiadiabatic 
regime u > A we find P = tanh[fi(w' - u)/(2k B T)}. 
This result is at first sight intriguing: for fast driving, 
the spin reaches thermal equilibrium like a static spin 
Hamiltonian characterized by an effective Larmor fre- 
quency (a/ — ui) = A[A/(2cj) + 0(A/uj) 3 } that vanishes 
for large uj. The polarization, which remains identically 
unity for all large frequencies at T = 0, decays eventu- 
ally at any finite T for large uj. Faster and faster driving 
at finite temperature enhances the pumped current up 
to u ~ HA 2 /ksT. For larger driving frequencies, ther- 
mal fluctuations catch up and suppress P causing the 



pumped current I to drop, as seen on the high-frequency 
side of Fig. [2j The reason why (u/ — uj) determines the 
Boltzmann occupancy of the two levels split by Huj' may 
be traced to the — uj ■ M term to be included in the ther- 
modynamically relevant functions for a body rotating at 
frequency uj, see e.g., §26 of Ref. [l^, where M is the 
body angular momentum which, for our spin, coincides 
with her. 

The results just presented are analytical and thus quite 
elegant and predictive. Obtained as they were for a 
rather special case however, how general are they? To 
address this question we solve numerically Eq. §5§ by 
means of Runge-Kutta integration [l8j |. and obtain ps{t) 
and from that I(t) = IoTr[cr y ps(t)]. The numerical ap- 
proach allows us to study the effect of unequal environ- 
ments in the x and z directions, where the simplifications 
leading to Eq. © do not hold. In particular, we con- 
sider the case a x ^ a z (by symmetry, it does not matter 
which one is larger). We find that at finite (but small) 
cvj the solution is no longer stationary even in the ro- 
tating reference frame chosen, and small oscillations of 
the density matrix and of the current at frequency 2uj re- 
main undamped in the long-time limit. Nevertheless for 
— > 0, the amplitude of these oscillating density-matrix 
terms vanishes linearly with , and the constant part of 
the density matrix at low temperature converges to the 
symmetric-environment case. In particular, at T = the 
polarization again saturates to 1. 

This behavior for — > can alternatively be re- 
covered by applying a rotating wave approximation to 
Eq. ([5]), i.e., by neglecting all the terms oscillating with 
frequency uj or uj' . Remarkably, the resulting equation 
again coincides without approximations with the one ap- 
propriate to the symmetric environment. We conclude 
that the results obtained for the symmetric environment 
are indeed representative of those expected in the more 
general asymmetric coupling case, provided the limit of 
weak coupling to the environment holds. In particular, 
Eq. ^ remains valid. 

The numerical solution of the master equation ([5]) also 
illustrates the transient approach to the stationary state. 
Figure [3] shows the full time evolution of the current, 
compared to the pure quantum evolution in the absence 
of dissipation. Note that, for a given coupling a, tem- 
perature affects not only the final steady current, via the 
final value of P, see Eq. (|6|) and inset of Fig. |3l but also 
the relaxation time with which this steady state is ap- 
proached in the initial transient. 

In Ref. 0] an investigation was attempted of nonadia- 
baticity, with numerical evidence that a stronger dissipa- 
tion might somehow compensate for the weak-coupling 
non-adiabatic current reduction relative to the geometric 
value of the adiabatic limit. Our exact solution clarifies 
that nonadiabaticity is fundamentally associated to such 
a radical current suppression that eventually, for large 
frequency, the charge pumped in one period Q p oc uj~ 2 . 
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a = 0.005. k B T/A = 0.01 
a = 0.005. k n T/A= 1 




FIG. 3: Time evolution of the current I(t) = I Tr[a v ps(t)] 
for cj = 2 A, in the dissipationless case (a = 0, solid curve), 
and in the transient induced by weak dissipation (a = 0.005), 
at low (dashed) and intermediate temperature (dot-dashed), 
starting from the initial \z, — } state. Inset: temperature de- 
pendence of the steady-state dc current /. 



Ref. [lg . 

In summary, we presented an analytical solution for 
the time-dependent pumping of DC current in a quan- 
tum model with dissipation, valid in the weak dissipation 
limit. The solution fully covers the crossover from the 
well-known adiabatic limit to arbitrarily high frequen- 
cies. The main physical surprise is that the frequency 
dependence of current is nonmonotonic, with an optimal 
value that moves from A upwards to infinity as tempera- 
ture is reduced. This effect, on the whole reminiscent of 
magnetic-resonance physics, could be directly detectable 
for example in multi-dot arrangements. 

We acknowledge helpful discussions with V. Brosco, 
R. Fazio, J. P. Pekola, and A. Russomanno. Research 
was supported by the Italian CNR through ESF Euro- 
core/FANAS/AFRI, by the Italian Ministry of University 
and Research, through PRIN/COFIN 20087NX9Y7, and 
by the French ANR, contract QNM ANR10-BLAN-0404- 
03. 



The method introduced in a very recent work [5J also 
can deal with the non-adiabatic regime (indeed our |n±) 
states coincide with the Floquet states of that method), 
and represents a more general approach to pumping prob- 
lems in the presence of weak dissipation 19]. However 
in that work intermediate frequencies were studied only 
in the absence of dissipation: the present model seems 
unique in affording an explicit analytic expression for the 
density matrix and current in the presence of dissipation 
for arbitrarily high frequency. 

Feasibility and quantitative estimates - Triple quan- 
tum dot systems have been recently realized experimen- 
tally [2(j, and could be used to implement the pump- 
ing effect proposed. By adopting a tentative hopping 
7o ~ 0.05 meV between the dots, we find a maximum 
current of the order of I — 0.05 meV q e /(V3h) ~ 7 nA. 
The magnetic field generated by this current could be 
detectable if the dot-ring arrangement had at least an 



[1] 



effective radius r c g ~ 200 nm [2l|. In that case, a ring- 
shaped SQUID of 5 /im radius placed ~ 5 above the 
quantum dots would intercept a flux of order 0.02 flux 
quanta, a value routinely detectable. 

The frequency and temperature regions where this cur- 
rent enhancement could be detected are determined by 
the frequency scale A of the effective spin-1/2 model. As- 
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be relevant in a molecular context as well, where an elec- 
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ample due to cyclic molecular distortions, see Fig. [T] and 
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